Quasistatic Scale-free Networks 
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A network is formed using the N sites of an one-dimensional lattice in the shape of a ring as nodes 
and each node with the initial degree k in = 2. N links are then introduced to this network, each link 
starts from a distinct node, the other end being connected to any other node with degree k randomly 
selected with an attachment probability proportional to k a . Tuning the control parameter a we 
observe a transition where the average degree of the largest node (fc m (a, N)) changes its variation 
from N° to iV at a specific transition point of ot c - The network is scale-free i.e., the nodal degree 
distribution has a power law decay for a > a c . 

PACS numbers: 05.40.-a, 64.60.Cn, 89.75.Hc, 89.75.-k 
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The nodal degree distribution function of a scale-free 
network (SFN) has a power law tail JjJ. Empirical data 
obtained for several social, biological and computational 
networks have confirmed existence of such probability 
distributions decaying as power laws j|, [| f|. For ex- 
ample the World-wide web S which is a network of web 
pages and the hyper-links among various pages and the 
Internet network || of routers or autonomous systems 
follow power laws as: P{k) ~ k^ 1 . The exponent 7 
varies between 2 and 3 for these networks. 

Networks are classified here as 'growing', 'quasistatic' 
and 'static'. In a growing network (GN) the nodes are 
introduced in the network one after another. After in- 
troducing the node a link is introduced connecting this 
node to one of its previous nodes. Therefore, in a GN 
both the numbers of nodes as well as links grow with 
time. In contrast, in the case of a quasistatic network 
(QN) a fixed number N of nodes are present at the ini- 
tial stage. N Links are then introduced one after another 
between pairs of nodes using some specific probability 
distribution. Therefore in a QN, the number of nodes in 
the network is fixed but the number of links grow with 
time. Finally, in a static network (SN) both the number 
of nodes as well as the number of links remain fixed and 
donot grow with time. 

Barabasi and Albert (BA) proposed Q a simple model 
for a growing SFN that has the following two essential in- 
gredients, namely: (i) A network grows from an initial set 
of m nodes with m < m links among them. Further, 
at every time step a new node is introduced and is ran- 
domly connected to m previous nodes, (ii) Any of these 
m links of the new node introduced at time t connects 
a previous node i with an attachment probability 7T,(t) 
which is linearly proportional to the degree fcj(t) of the 
i-th node at time t: irf A (t) ~ k. t (t). For BA model 7 = 3 

B- 

The second criterion reflects the phenomenon of 'rich 
gets richer' i.e., a node with a large degree attracts more 
nodes to get linked. Krapivsky et. al. showed that this 
linear dependence is a necessary condition and any other 
non-linear dependence on the degree destroys the scale- 




FIG. 1: A quasistatic scale-free network with 128 nodes 
placed on a circular ring so that the initial deg ree k{ n n of 
each node is equal to 2. Nodes are then selected one after 
another without repetition and are randomly linked to other 
nodes with an attachment probability proportional to k a . In 
this picture a = 2.6 is used and a large degree node is visible. 



free behaviour i.e., the power law variation of the degree 
distribution 0. Indeed, if in general the degree depen- 
dence has a variation with the a-th power of the degree, 
it has been shown that the scale-free nature of the BA 
model network exists only for a — 1 and for no other 
value of it, smaller or larger. 

Several other interesting networks are also studied. 
Erdos and Renyi studied long ago random graphs of 
N nodes where links are introduced with probability 
p between arbitrary pairs of randomly selected nodes 
M. A largest component of such a network connecting 
nodes of the order of N appears at a particular value of 
p c = 1/N. Very recently SFNs are studied on Euclidean 
spaces where the BA attachment probability is modified 
by a link length I dependent factor || [l(| O, Load 
distribution in Scale-free networks has also been studied 
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FIG. 2: The variation of the order parameter <f>(a, N) which 
is the average degree of the largest node divided by N with a 
for three different network sizes N = 2 8 (circle), 2 10 (square) 
and 2 12 (triangle). The inset shows the plot of the threshold 
values a c (N) with N~ ' 2 for seven different TV values from 2 7 
to 2 13 increasing by a factor of 2 at each step. 
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In this paper we ask the question if the growing con- 
dition of the BA model is really a necessity to achieve 
a scale- free network. We will see below that it is not, a 
suitable choice of the attachment probability in a qua- 
sistatic network may result a power law decay for the 
degree distribution as well. We call our model as the Qu- 
asitatic Scale- free Network (QSFN). Recently Doye has 
shown that the network topology of a potential energy 
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landscape is a static scale-free network |14 



a quenched intrinsic fitness to every node and using the 
attachment probabilities depending on the fitnesses it is 
also possible to get SFNs without growth and preferen- 
tial attachments |1| . A steady state model for scale- free 
graphs has also been proposed Jig . 

Our quasistatic model starts with an one-dimensional 
regular lattice in the form of a circular ring as in the 
Watts and Strogatz's Small- world Network (n). The lat- 
tice has N nodes serially marked from i—1 to N and 
N links between successive pairs of nodes with periodic 
boundary condition. Each node is therefore connected to 
only two of its nearest neighbours situated on the oppo- 
site sides and the degree ki for each node i is exactly 2 
to begin with. We add to this system another N distinct 
links in total, such that a new link starts from each node. 
A t-th. link is added at time t and one end of it is attached 
to the t-th. node, the other end of the link is connected 
to a node j with degree kj(t) selected randomly from the 
rest of the N — 1 nodes of the system using the following 



FIG. 3: The local derivative of the order parameter 
d<f>(a, N)/da has been plotted with a for 7 different N values 
starting from 2 7 and increasing by a factor of 2. Scaling of 
the derivative is done in the inset for N = 2 8 , 2 10 and 2 12 by 



plotting [d<j)(a,N)/da]N~ 
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attachment probability: 



(1) 



where a is a continuously tunable parameter. Therefore 
when the network is complete no node is left with degree 
2. 

At this point we would like to clearly distinguish be- 
tween our QSFN and the model B Q. Model B also 
starts with a collection of N nodes but no links, so that 
the initial degree of each node is strictly zero and the 
graph has N components each with only one node. At 
each time step a node i is selected randomly with uni- 
form probability and is connected to another node j with 
a probability TVj(t) ~ kj(t). This implies that isolated 
nodes are being linked one by one to a single connected 
component of the graph. For this component, both the 
number of nodes as well as the links grow with time and 
therefore model B is clearly a growing network. This is in 
contrast to our QSFN model where the initial network is 
a N node connected graph which grows by adding further 
links but not the nodes any more. 

In the case of a = in QSFN, all nodes have equal 
probabilities to get connected by a new link and this is 
the random network model || for which the degree dis- 
tribution is the Poisson distribution [0. For a > 0, a 
node with larger degree has a higher probability to get 
connected by links. The case of a = 1 corresponds to 
the similar attachment probability as in the BA model of 
SFN. Our simulation results show that the degree distri- 
bution in this case has an exponentially decaying form: 
P(k) ~ exp(— ak) with a « 1.08. On continuously in- 
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FIG. 4: Scaling of the order parameter: [4>{oi, N) — 
<P(1,N)]N - 025 is plotted with [q - a c (N)]N 1/3 for the data 
collapse. The inset shows that average degree of the largest 
node (fe m (a, N)} varies as iV M at the transition point a c where 
we obtained fi — 0.91 ± 0.04. The dotted line, having a slope 
unity, is given to compare with the data. 



creasing a further larger degree nodes become more prob- 
able and we see that the degree distribution decays less 
sharply and changes to a stretched exponential form as: 
P{k) ~ exp[— fc x ( Q ))]. The exponent x( a ) decreases con- 
tinuously from its value 1 at a — 1 to \ ~ 0.3 at a = 1.75. 

On increasing the parameter a further the system 
makes a transition to a different behaviour where a single 
node having the maximum degree k m connects to a finite 
fraction of the N nodes. However this transition takes 
place at a specific value a c of a so that the degree distri- 
bution has a power law tail and therefore the network is 
scale- free for all a > a c . 

Naturally the control parameter in this model is a 
where as like the phenomenon of percolation, an order 
parameter in this problem may be the average degree of 
the largest node <f)(a,N) = (k m (a,N))/N. We monitor 
the variation of (f>(a, N) with a and plot them in Fig. 
2 for networks of three different sizes TV = 2 8 ,2 10 and 
2 12 . A sharp increase in the order parameter is observed 
around a = 2, again the sharpness of the curves increases 
with TV. The first derivative d(p(a, N) / da of the order pa- 
rameter is plotted in Fig. 3 for seven different network 
sizes. Every curve has a peak at some N dependent value 
of a c (N) where <j>(a, N) increases at the fastest rate. For 
locating the precise value of a c (oo) where this transi- 
tion is taking place for an infinitely large network, we 
plot a c (N) values as a function of N^ 1 ^ in the inset of 
Fig. 2. Using a trial value of v — 5 we could get all 
seven points on a straight line and on extrapolating to 
the N — > oo limit we get a c (oo) = 1.85 ± 0.10. 



All the first derivative curves are suitably scaled in the 
inset of Fig. 3 and the data are collapsed. Each curve is 
first shifted by an amount a c (N). Then the abscissa is 
scaled by TV 018 and the ordinate is scaled by jV~ 0,21 and 
we get a nice collapse of three curves of sizes N = 2 8 , 2 10 
and 2 12 . This analysis implies that the order parameter 
curves in Fig. 2 become progressively sharper as the 
network size N gradually increases and 4>(a, N) increases 
very rapidly at a c (oo). This variation is very similar to 
the variation of the order parameter (the fraction of mass 
in the largest cluster) in the percolation phenomena. 

The scaling of the order parameter <p(a, N) is shown 
in Fig. 4. Asymptotically as N — > oo the <fi(a, N) — > 
for a < a c but for finite N we subtract the </>(l, N) from 
4>(a, N) and try to scale only the difference. We assume 
the following scaling behaviour: 

[<t>(a, N) - 0(1, N)]NM" = T[(a - a^N 1 ^} (2) 

where the scaling function T(x) — > x 13 for x « 1. From 
Fig. 4 we get v = 3 and (3 — 0.075. This value of v is 
not very consistent with its value 5 obtained in the inset 
of Fig. 2. We believe this difference is due to finite size 
of our simulations. 

The dependence of the average degree (k m (a c , N)} of 
the largest node right at the transition point a = a c is 
also studied and plotted in the inset of Fig. 4 on a double 
logarithmic scale. Assuming a power law dependence on 
the network size with an exponent /i as (k m (a c , N)) ~ 
the plot gives a value of /i = 0.91 ± 0.04. Here /j, 
is an exponent similar to the fractal dimension of the 
infinite incipient percolation clusters at the percolation 
threshold. 

Finally the degree distribution P(k,a c (N),N) is cal- 
culated for different network sizes N at their transition 
points a c (N). In these calculations we have not counted 
the maximum degree node similar to what is usually done 
to estimate the cluster size distribution at the percolation 
threshold. Except for very small k values all distribu- 
tion curves have straight portions on double logarithmic 
plots and the region over which nearly constant slopes 
are observed, increases with increasing N. The fluctuat- 
ing data is suitably binned: for each k value the data is 
averaged over the bin from k to 2k — 1 and is plotted at 
[k(2k- l)] 1 / 2 . We assume a scaling form like: 

P(fc,a c (A^),A)A' ) ~g(k/N*) (3) 

where Q (x) is expected to be an universal scaling function 
such that Q(x) — ► x 1 for x << 1 so that 7 = 77/^. The 
scaling plot is shown in Fig. 5 where we get 77 = 3 and 
£ = 1 giving 7 = 3 as in the BA model of SFN. This is a 
very interesting result that even on a static network with 
a suitable attachment probability k a with a non-trivial 
value of a one gets a scale-free degree distribution for 
the network, the distribution being the same as the BA 
model of SFN. 
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ing to each node. That means we systematically attach 
the one end of each link to one node, the other end of 
the link is probabilistically attached to any other node 
of degree k with a probability proportional to k a . Our 
numerical study indicates that there exists a transition 
point a c beyond which the resulting network has a scale- 
free structure so that the degree distribution has a power 
law tail for a > a c . 
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FIG. 5: Scaled degree distributions at the transition points 
a c (N) for different network sizes have been plotted. Data 
collapse is obtained for large k values when P(k, a c , N)N V 
has been plotted with k/N^ with r/ — 3 and £ = 1. This gives 
the degree distribution exponent 7 = r//^ = 3. 



What happens when a > a c l Our observation is that 
degree distribution is still a power law but the exponent 
7 is now a dependent. In the range of a > 3, the degree 
exponent 7 is approximately equal to a whereas for a c < 
a < 3, 7(a) ~ 3. 

We also have tried a more stochastic version of this 
model where to add a link, we select in parallel two nodes 
arbitrarily from the whole set of N nodes using the at- 
tachment probability in Eqn. (1). If these two nodes are 
not the same node, we connect them. This is also a QSFN 
but the difference here is both the nodes of each link are 
selected randomly where as in the QSFN described above 
only, one node is selected randomly and the other node is 
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sele cted systematically one after another from the sot of 



nodes. Our numerical studies indicate that behaviours of 
both versions are very similar. 

Another variation of our model has been studied with 
initial degrees of all nodes as fcj n = 1 which corre- 
sponds to a situation where alternate pairs of nodes on 
an one-dimensional ring shaped lattice are linked. Again 
we numerically observe that the a c — 1.67 ± 0.10 and 
■y(a c ) w 2.5 for this model. 

To summarize, we have studied a quasistatic scale-free 
network (QSFN) where we have a collection of TV nodes 
initially present, each node is having k in initial degree. A 
system of TV links are then introduced, one correspond- 
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